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As a possible way of modifying the intrinsic properties of graphene we study the doping of graphene 
by embedded boron clusters with density functional theory. Cluster doping is technologically relevant 
as the cluster implantation technique can be readily applied to graphene. We find that Br clusters 
embedded into graphene and graphene nanoribbons are structurally stable and locally metallize 
the system. This is done both by the reduction of the Fermi energy and by the introduction 
of boron states near the Fermi level. A linear chain of boron clusters forms a metallic "wire" 
inside the graphene matrix. In a zigzag edge graphene nanoribbon the cluster-related states tend 
to hybridize with the edge and bulk states. The magnetism in boron doped graphene systems is 
generally very weak. The presence of boron clusters weakens the edge magnetism in zigzag edge 
graphene nanoribbon, rather than making the system appropriate for spintronics. Thus the doping 
of graphene with the cluster implantation technique might be a viable technique to locally metallize 
graphene without destroying its attractive bulk properties. 

PACS numbers: 73.22.Pr,75.75.-c,71.55.-i,75.75.Lf 



I. INTRODUCTION 

The two-dimensional (2D) honeycomb lattice of 
graphene is the basic building block of OD fullerenes, ID 
nanotubes or 3D graphite. Besides being a 2D topologi- 
cal structure with high strength and flexibility, graphene 
exhibits fascinating, and rather unique electronic prop- 
erties, such as behaving like massless Dirac fermions and 
showing a particular high charge carrier mobility and 
electrical conductivityi^. These exceptional properties 
lead to huge scientific interest in graphene, and they make 
this particular material a promising candidate for future 
technological applications^"—. 

A free-standing graphene sheet is a semimetal. How- 
ever many applications require semi-conducting proper- 
ties. A band gap can be induced in graphene systems, 
for example by oxidation or by considering graphene 
nanoribbons (GNRs). Narrow width ribbons (< 10 nm) 
have a band gap due to lateral quantum confinement 
and edge effect s^^^^. A process for fabricating dense 
graphene nanoribbon arrays with widths as narrow as 
10 nm has been developed quite recently using self- 
assembled patterns of block copolymers on graphene, 
which was grown epitaxially on SiC. 

The edge geometries of GNRs can be rather diversei^. 
Usually armchair or zigzag types (AGNR or ZGNR) 
are considered in the literaturei^ — . AGNRs are ei- 
ther metallic or semi-conducting depending on the rib- 
bon width, and ZGNRs have metallic properties with 
flat energy bands near the Fermi level Ep. The lat- 
ter give rise to a sharp peak in the electronic density 
of states (DOS) right at Ej?. It is currently believed that 



this peak at is avoided by magnetic interactions, and 
that the ground state of ZGNRs is actually an antifer- 
romagnetic semi-conductor (AFM, atoms are ferromag- 
netically ordered along one edge and antiferromagneti- 
cally ordered between opposite edges) with a small band 
gap. In density functional theory (DFT) calculations the 
AFM magnetic ground state is almost degenerate to a 
ferromagnetic (FM) stat o^^i^^ . These results suggested 
the use of ZGNRs in spintronic applications, and they in- 
spired tremendous theoretical work. It is understood now 
that there should either be some translational symmetry 
breaking between edges in GNRs or some local sublat- 
tice imbalances in graphene/graphite to create localized 
states at Ep to obtain a finite magnetic moment. Several 
approaches have been proposed, such as the introduction 
of vacancies^^"— , networks of point defects^, impuri- 
ties/chemical dopin g^^i^^ — , partial hydrogenation^*^, 
edge roughness^^, electron or hole doping^ or an ap- 
plied electric field to introduce half- metallicity^^. How- 
ever, some recent studies also point out that the edge 
states and the associated edge magnetism might be ei- 
ther very weak, or absent, or even artefacts of improper 
theoretical methodolog}*^ — . 

Another way of modulating and controlling the elec- 
tronic properties of graphene is chemical doping, which 
was widely studied in the pasi,29|32 42^^^ rj.^^ 

synthesis of substitutionally doped graphene samples is 
presently reported in the literature, such as B- and N- 
doping by arc discharge^^, N doping by chemical vapor 
deposition (CVD)^ methods, and doping with transi- 
tion metals atoms^. B/N^^, B^^^, P^, transi- 
tion metals atoms^ doped graphene systems were also 
investigated from a theoretical point of view, with par- 
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ticular focus on the effects of substitutional doping on the 
structural and electronic properties of graphene. Boron, 
the neighbor of carbon in the periodic table, seems to 
be an ideal dopant^^. However, individual boron atoms 
in graphene are rather mobile and tend to aggregate at 
the edges and not in the bulk of the sample, as stud- 
ies of boron doped graphene nanoribbons^^ and boron 
doped carbon nanotubes^^ show. In the present article 
we do not suggest to dope graphene with individual boron 
atoms, but with entire boron clusters. 

Small boron clusters have a quasiplanar geometry and 
can be viewed as an assembly of pyramidal build- 
ing blocks (Aufbau principle}^. Boron and carbon have 
roughly the same covalent radius. Therefore, one may 
expect that boron could easily be integrated into carbon 
systems. On the other hand, nitrogen adsorption causes 
some structural distortions on graphene, despite having 
a similar atomic radius to carbon as well^^. 

We have previously reported that boron cluster can 
be easily accepted by graphene^^. Ion beam meth- 
ods can be utilized efficiently to dope boron into the 
graphene/graphite matrix, either as atomic site a^^i^^ or 
as cluster^. A cluster implanter particularly suited for 
boron (5io, Bis clusters) has been developed with the en- 
hancements of ultra low energy beam currents^. Studies 
of doping of few-layered graphene and carbon nanotubes 
for a range of dopants (including boron) using ion im- 
plantation techniques^ , as well as atomistic simulation 
of low energy boron and nitrogen ions implantation into 
graphene^^ have recently been reported in the literature. 

In the present work we study the effect of embedded 
boron clusters on the electronic and magnetic properties 
of graphene systems. We consider 2D graphene and ID 
AGNRs and ZGNRs. These systems are substitution- 
ally doped with B7 clusters to create (1) localized states, 
(2) local sublattice imbalances and (3) localized magnetic 
moments. We find that the B7 clusters locally metal- 
lize the graphene matrix, and that linear chains of boron 
clusters form metallic wires inside the graphene matrix. 
Furthermore the presence of the B7 clusters in magnetic 
ZGNR weakens the magnetism in such systems. 



II. COMPUTATIONAL METHOD 

Our calculations are based on the density func- 
tional theory^ (DFT) within the framework of the 
generalized gradient approximation (GGA) using the 
Perdew-Wang parametrization (PW91)^i^. We used 
the Vienna ab initio simulation package (VASP, ver- 
sion 4.6)^1^ employing the projected augmented wave 
(PAW) method^UiTi. 

In our study, we consider graphene super cells, 9- 
AGNR(+H) and 10-ZGNR(+H)s and their boron doped 
counterparts (+BD, +isolated-BD). The prefix "9 (10)" 
represents the number of armchair (zigzag) lines cross the 
width of the ribbon, and the suffix "+II" indicates that 
each carbon edge atom is passivated by one hydrogen 



atom. A suffix "+BD" that the structure is boron doped 
and "+isolated-BD" indicates that the embedded boron 
clusters in adjacent unit cells of the graphene supercell 
systems are separated by at least 10 A. 

A supercell system still repeats itself periodically in 
both the horizontal (x-) and the vertical (y-) direction. 
In order to simulate GNRs, and in order to avoid any 
residual interactions, each replica is separated from its 
neighboring replica by least 10 A in both, the edge-edge 
(to the right or to the left) and the (z-) direction perpen- 
dicular to the planes of the replicas^^. 

A plane wave basis set with a kinetic energy cutoff 
value of 450 eV is used for all calculations. For the k- 
point sampling of the Brillouin zone, a F-point centered 
Monkhorst-Pack grid was used. The k-space integration 
was carried out with the method of Methfessel and Pax- 
ton in first order and a smearing width of 0.1 eV. Op- 
timal k-point meshes are individually converged for each 
system by reducing the changes in the total energy at 
least below 10 meV/atom^^ After obtaining the opti- 
mal k-point meshes, the geometries of all systems were 
fully optimized (relaxed) in a self-consistent fashion. The 
relaxed structures were close to the unrelaxed ideal ge- 
ometries. After the relaxation of each structure an addi- 
tional static self-consistent calculation was performed to 
obtain the charge density and the total energies of the 
system, where /c-space integrations were carried out us- 
ing the tetrahedron method^^. In all of our self-consistent 
calculations the total energies were converged, such that 
energetic changes were less than 10 meV. Band structures 
and DOS of the optimized systems were computed non- 
self-consistently by using the charge densities obtained 
in the final self-consistent calculation. At least twice the 
optimal size of /c-point meshes was used for the DOS cal- 
culations. For magnetic systems collinear, spin polarized 
(SP) DFT calculations were performed in all the steps 
described above. The different magnetic states were de- 
termined by providing different initial spin densities at 
the beginning of the self-consistency cycles. System sizes, 
total energies, and magnetic moments of the systems con- 
sidered for this study are summarized in table HI In the 
following, the non-spin-polarized calculations will be ab- 
breviated by NSP. 



III. RESULTS AND DISCUSSION 

In the presented study, a seven atom planar boron B7 
cluster is embedded into graphene and various GNRs (see 
for example Fig. HJ^h)). The basic stability of this struc- 
ture was reported by us in an earlier stud y^^i^^ . The 
boron atom, being one electron short in comparison to 
the carbon atom, is a p-type dopant. Embedding a B7 
unit (21 valence electrons) into a graphene system leads 
to the replacement of a Cq hexagon (24 valence elec- 
trons), which results in a decrease of 3 electrons (intro- 
duction of 3 holes) per unit cell, and thus in the lowering 
of the Fermi energy Ep. The degree of doping is calcu- 
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lated with the fohowing equation 

A^F = Ef'"'' - (1) 

and numerical values for the shift of the Fermi energy are 
tabulated in Table ID 

The atomic structures together with partial charge 
densities (PCD), projected DOSs, plain band structures 
and spin densities of an optimized graphene supercell and 
of GNR+H systems with boron cluster doping are shown 
in Figs. [T]|4l^a)-(d), respectively. The PCD is defined as 

PCD(r) = ^ r^^^E 5{E - Eni.Mni.{v)\\ (2) 

where 0nk(r) are the Kohn-Sham orbitals with the band 
index n and the wave vector k, Env. are the band ener- 
gies, the summation over k represents a Brillouine zone 
average, and the parameter e = 0.1 eV in all our calcula- 
tions. The PCD visualizes the localization of electronic 
states and it is drawn on top of the atomic structure. 
Note that the chosen energy range is relevant for electric 
conduction. 



A. Localization of Metallicity 

1. Isolated Boron Cluster in Graphene 

First we study the properties of an nearly isolated B7 
cluster in graphene (graphene+isolated-BD system). We 
use a relatively large supercell to ensure that the periodic 
replica of the clusters do not interacti^. The optimized 
structure together with the PCD in the specified energy 
window is depicted in Figure [IJa). The PCD is basi- 
cally localized on the outer boron atoms of the cluster 
and on nearby carbon atoms, whereas the central atom 
does not contribute to the states near Ep- These elec- 
tronic states have tt character, as seen from the projected 
DOS plot in Figure [TJc). The a states are responsible 
for in-plane stability of boron clusters within a graphene 
matrix, and they do not contribute to the states near 
Ep' The general form of the DOS of the pristine and the 
boron doped systems in Figs. mb,c) is quite similar, in- 
dicating that the boron clusters do not significantly alter 
the bulk properties of the system in this energy window. 
Essential changes happen near Ep^ where the DOS has a 
pronounced peak. The peaks indicates that the induced 
boron states are localize, which is also confirmed by look- 
ing at the PCD. 

The change of the band structures from Fig. [Hd) to 
(f) can be described as twofold: (i) creation of new states 
near Ep that are primarily localized on the boron atoms 
(as is obvious from the PCD plot shown above) and (ii) 
reduction of Ep by A£^f ca. 0.2 eV (see Tab. |T]) in a 
way, that cannot be described by the rigid band model. 
Though the system has a finite DOS(£^f), and despite 
the fact that we see bands crossing Ep in Fig. ffl^f), the 
system is not necessarily a metal, because the presumed 




FIG. 1. (Color online) Electronic structure of graphene with 
embedded B7 clusters that are (nearly) isolated from each 
other (graphene+isolated-BD); the center-to-center separa- 
tion of the boron clusters in the y-direction is 17 A. (a) Non- 
spin-polarized (NSP) partial charge density (PCD) for states 
in the energy window \E — Ef\ < 0.1 eV at contours of 0.01 
e/A^; the PCD is basically localized on the outer boron atoms 
of the cluster as well as on nearby carbon atoms. (b,c) NSP 
density of states (DOS) of pristine and boron doped systems, 
respectively; orange (light gray) represents the tt character 
and blue (dark gray) represents the a character, (d-g) Band 
structures for NSP pristine, spin-polarized (SP) pristine, NSP 
boron doped and SP boron doped cases, respectively, (h) The 
spin density at contours of 0.01 /xb/A^ is localized on the clus- 
ter. In (a) and (h) four unit cells are depicted. 
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conduction channels near Ef are highly localized and 
have little overlap. However we can conclude that the 
boron cluster locally metallizes the system. In order to 
create pronounced metallic conduction channels one has 
to decrease the distance between the clusters. This is 
what we will do in the next subsection. 



2. Boron Cluster Chain in Graphene 

In the graphene+BD system the distance between the 
clusters along the y-direction is halved in comparison to 
the previous case. The atomic structure along y can be 
described as an alternating By-Ce chain. The geometry 
and the corresponding PCD are depicted in Figure [2fa). 
The PCD indicates that the states near Ef form a ID 
metallic "wire" inside the 2D graphene matrix. The over- 
lap between states related to the cluster tends to increase 
the dispersion of the the bands near Ep in Fig. [2jf)). 
This overlap also reduces the peak in the DOS{Ef) in 
Fig. [2fc). As the cluster density is doubled compared 
to the previously considered case the doping level is also 
approximately doubled to AE-p ca. 0.4 eV. The presence 
of boron a states in the occupied part of the DOS (at -1 
and -3 eV) is qualitatively similar to the previous case, 
but this time it is much more pronounced. 



3. Boron cluster chain in AGNR 

The PCD of the alternating By-Ce chain in an 9- 
AGNR+H is shown for the optimized geometry in Figure 
[Sja). A comparison with Figure [2]^a) indicates that the 
PCDs of the two systems are almost identical. Thus the 
finite size of the AGNR does not influence the basic elec- 
tronic properties of the system close to Ep- This further 
supports our finding that the metallization of the con- 
sidered graphene systems by the boron clusters tends to 
happens on a very small length scale. The undoped 9- 
AGNR+H system is a semi-conductor with a band gap of 
0.37 eV (see Fig.j^d)). Doping and appearance of boron 
related tt and a states turn the system into a ID metal 
with 4 (3) bands at Ep for the NSP (SP) cases (this band 
counting includes the spin degeneracy). The a states at 
about -1 eV in Fig.[3fc) are now much more smeared out 
in energy compared to the previous systems. 



4. Boron cluster chain in ZGNR 

The undoped 10-ZGNR+H exhibits a considerable 
peak in the NSP DOS(^f) (see Fig. H^b)) coming from 
the flat bands at Ep in the NSP band structure (see 
Fig. m^d)). This peak is associated with the edge states 
mentioned in the Introduction. As discernible in Fig- 
ure IHa), introducing the alternating By-Ce chain in a 
10-ZGNR+H system leads to a dramatically different 
electronic structure in comparison to the previous cases. 




(h) 



I I 




FIG. 2. (Color online) Electronic structure of graphene with 
B7 clusters forming a linear chain with a center-to-center sep- 
aration of 8.5 A(graphene+BD). (a) The NSP PCD indicates 
that the states near Ep form a ID metallic "wire" inside the 
2D graphene matrix. (b,c) NSP DOSs of pristine and boron 
doped cases, respectively, (d-g) plain band structures of NSP 
pristine, SP pristine, NSP boron doped and SP boron doped 
cases, respectively, (h) The spin density shows that the mag- 
netic moments on the different clusters are smaller than the 
ones on the isolated clusters in Fig. [T] Eight unit cells are 
depicted in (a) and (h). Color codes, energy ranges, contours 
and abbreviations are the same as in Fig. [T] 



There are three qualitative types of states that contribute 
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FIG. 3. (Color online) Electronic structure of B7 clusters 
forming a linear chain within an armchair graphene nanorib- 
bon that has hydrogen terminated edges (9-AGNR+BD+H) 
(a) The NSP PCD is almost identical to the one in Fig.[2la), 
i.e., the states form a metallic "wire" inside the nanoribbon. 
(b,c) NSP DOSs of pristine and boron doped cases, (d-g) 
Plain band structures for NSP pristine, SP pristine, NSP 
boron doped and SP boron doped cases, respectively, (h) 
spin density at contours of iso = -0.01 /xb/A^ Three unit 
cells are depicted in (a) and (h). Color codes, energy ranges, 
contours and abbreviations are the same as in Fig. [TJ 



to the PCD: localized boron states, delocalized carbon 
bulk states and localized carbon edge states. The pres- 
ence of the boron clusters allows both the carbon bulk 
states and the edge states to hybridize with the cluster 
induced states. This results in the rather complex and 
very delocalized PCD shown in Fig. HJ^a) 




FIG. 4. (Color online) Electronic structure of B7 clusters 
forming a linear chain within a zigzag graphene nanoribbon 
that has hydrogen terminated edges (10-ZGNR+BD+H) (a) 
The NSP PCD shows that the states from the boron cluster 
hybridize with the edge and bulk states. (b,c) NSP DOSs 
of pristine and boron doped cases, (d-g) Band structures for 
NSP pristine, antiferromagnetically (AFM) SP pristine, NSP 
boron doped and ferromagnetically (FM) SP boron doped 
cases, respectively, (h) The spin density of the AFM state 
shows that magnetic moments on and in the vicinity of the 
cluster are suppressed. Thus the presence of the B7 clusters 
weakens the edge magnetism. Three unit cells are depicted 
in (a) and (h). Color codes, energy ranges, contours and 
abbreviations are the same as in Fig. [T] 



Due to the interplay of the three different types of 
states the change of the electronic structure in the DOS 
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and band plots in Fig. |4]is relatively complicated. But, 
similar to the previous cases, the DOS near Ef increases, 
and there are more bands near Ep than in the undoped 
case due to the appearance of the boron related bands. 
In the AFM state the boron doping increases the num- 
ber of bands at Ep from to 4, and in the FM state 
the increase is from 2 to 4. The a states near -1 eV in 
the DOS of the doped system in Fig. HJ^c) are most dis- 
persed in the present system, and they come very close 
to Ep- Further doping of the system could eventually 
create holes in these kinds of a bands. It is interesting 
to note that the superconductivity in MgB2 and boron 
doped diamond is related to holes in a bands^^^^^. One 
might speculate whether boron doped ZGRNs could al- 
low for a similar type of superconductivity. 

To conclude this section, we found that By clusters 
embedded into graphene systems are structurally stable, 
and they strongly modify the electronic structure near 
Ep. The boron clusters locally metallize the system. 
This is achieved both by the reduction of Ef^ and by 
the introduction of boron states near Ep. A linear chain 
of boron clusters forms a ID metallic "wire" inside the 
graphene matrix. In a ZGNR the cluster-related states 
hybridize with the edge and bulk states. 

B. Energetics and magnetic properties 

In this section we will investigate the energetics and 
magnetic properties of the boron doped graphene sys- 
tems. The results comprising the energies and the mag- 
netic moments for nonmagnetic (NM) and magnetic (SP, 
FM, AFM) states are tabulated in Table II The ener- 
gies of different magnetic states are characterized by the 
difference in their total energy as compared to the total 
energy of the magnetic ground state (E^f): 

A£^mag = (^tot - E^o^)/NmA, (3) 

where Nma is the number of magnetic "centers" per unit 
cell. For the majority of systems under consideration 
A^MA = 1, except for the ZGNR systems, where this 
value is equal to the number of carbon edge atoms per 
unit cell (A'ma = 6). When analyzing the spin densi- 
ties in Figs. mh)-|ll^h) it becomes quite obvious that the 
magnetic moments are never concentrated on a single 
atom. So A'ma represents the number of similar groups 
of atoms per unit cell, which collectively carry a certain 
total magnetic moment. 

Undoped graphene and AGNRs are nonmagnetic. 
However, when graphene is doped with an isolated boron 
cluster (graphene+isolated-BD) the system has a total 
magnetic moment of 1 /i^ per cluster due to the odd 
number of electrons per boron cluster (21 in B7 as com- 
pared to 24 in Ce). Thus a boron cluster carries an 
unpaired electron that has an magnetic moment of ex- 
actly 1 fiB' The corresponding spin density shown in 
Fig.fljh) is primarily concentrated on the outer B atoms, 
and the latter does not interact with the cluster replica 



TABLE I. The energies, magnetic states and moments of all 
systems considered in this study. "State" denotes different 
systems and magnetic orderings: UC (unit cell), SC (super- 
cell), H (hydrogen terminated), BD (boron doped), NM (non- 
magnetic), SP (spin polarized), FM (ferromagnetic), AFM 
(antiferromagnetic). Aat is the number of atoms per unit 
cell, /i is the total magnetic moment per unit cell, ^tot is the 
total energy per unit cell. The magnetic ground state is indi- 
cated in boldface. A^^mag is defined in Eq. [3] and the degree 
of doping A^F is defined in Eq. [U 
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in the adjacent unit cells. By bringing the boron clus- 
ters closer to each other within the graphene-BD linear 
chain system, the magnetic moment is quenched down 
to 0.3 /i^/cluster, due to the onset of magnetic exchange 
interactions between the clusters. The reduction of the 
magnetic moment is also discernible in the spin density 
in Fig.[2fh). The magnetic properties do not change very 
significantly from the graphene-BD linear chain system 
to the 9-AGNR+BD(+H) systems. The magnetic mo- 
ment is still 0.3 /i^/cluster, and the spin densities of the 
two systems are relatively similar. This again indicates 
the the change of the electronic and magnetic structure 
induced by the presence of the cluster (s) is of a rather lo- 
cal nature. However there is one thing to be pointed out: 
the ground state of the 9-AGNR+BD system is magnetic 
and the ground state of the 9-AGNR+BD+H system is 
NM. We speculate that the hydrogen passivation removes 
dangling bonds at the edges, and it also minimizes the 
exchange interactions that let the system prefer to settle 
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down in a NM ground state. However, the energy differ- 
ences between the non-magnetic and the magnetic states 
(AEmag) in the graphene and AGNR systems are very 
smah, i.e. between 6-9 meV/atom (Table |T]). Such smah 
values indicate that the magnetism in these systems must 
be very weak. 

ZGNRs are believed to carry magnetic moments along 
the edges. However the magnetic interaction between the 
two edges of the ribbon is very weak, thus leading to a rel- 
atively unstable magnetic system^. As shown above, the 
boron cluster usually carries a magnetic moment. So one 
might speculate that the magnetic interaction between 
the edges in ZGNRs could be enhanced by the presence 
of the boron clusters. The more as both the cluster mag- 
netic moments and the ones on the edges are generated 
by TT-electrons, and therefore the two systems should be 
able to interact. 

The answer to this question is given in Fig. Hl^h). It 
shows that the edge states are weakened by the presence 
of the boron cluster. The boron cluster itself does no 
longer carry a net magnetic moment. Thus in ZGNR 
systems the boron cluster locally supresses the magnetic 
moments, and it effectively weakens the magnetic inter- 
actions. Note that the weakening of the edge magnetism 
happens both in the FM and in the AFM state. Fur- 
thermore, for the different FM states that are listed in 
Tab.Uthe total magnetic moment is always smaller in the 
doped case than in the undoped case. In the 10-ZGNR 
system the boron cluster is located at a very asymmet- 
ric position in order to break the symmetry between the 
edges. As result of the asymmetric supression of mag- 
netic moments the positive and negative moments of the 
AFM states on opposite sides of the ribbon are no longer 
equal in magnitude, and there is a residual total mag- 
netic moment per cell of 0.1 /iB/cell and 0.2 /i^/cell for 
the unpassivated and H passivated cases, respectively (as 
opposed to fiB /cell for the undoped systems). In the 
undoped systems the AFM and FM states are energeti- 
cally separated by only AEmag = 5 — 6 meV/atom. In 
the doped systems the two states are isoenergetic with 
A^mag = 1 meV/atom (smaller than the DFT energy 
resolution). This again is a strong indication that mag- 
netic interactions are suppressed by the presence of the 
boron cluster. Note that a similar suppression of the edge 
magnetic moments with increasing Ti concentration near 
the edge of ZGNRs is reported in Ref.^^. 

We can summarize this section by saying that the mag- 
netism in boron doped graphene systems is generally very 
weak. Furthermore the edge magnetism in ZGNRs in 
weakened by the presence of boron clusters. 

IV. SUMMARY AND CONCLUSION 

We studied the doping of graphene and various 
graphene nanoribbons by embedded B7 clusters. The 
embedding seems to be stable, and the boron clusters 



tend to strongly modify the electronic structure near Ep, 
which leads to a very local metallization of the graphene 
substrates. This metalliation is achieved by both the re- 
duction of Ef^ and by the introduction of boron states 
near Ep. 

We also showed that a linear chain of boron clusters 
forms an effective ID metallic "wire" inside a semicon- 
ducting substrate made of a suitable graphene nanorib- 
bon, where the metallicity rapidly decays away from the 
wire. Given the fact that a semiconducting graphene sub- 
strate does not require any additional doping, any con- 
trolled layout of linear chains of boron clusters might ac- 
tually be the technological basis of a layout of integrated 
circuits, which is at least one or two length scales below 
state-of-the-art microelectronics. Note that the corre- 
sponding cluster implantation technique can already be 
applied to graphene. 

Furthermore the presence of the boron clusters also 
introduces a states in the occupied part of the band 
structure within 1 eV below the Fermi level. We wonder 
whether this particular feature could give rise to electron- 
phonon mediated superconductivity, similar to MgB2 or 
boron doped diamond. 

Finally, the magnetism in boron doped graphene sys- 
tems is generally very weak, and the edge magnetism 
in ZGNRs is further weakened by the presence of boron 
clusters. Therefore we do not foresee any application of 
embedded boron clusters in the field of graphene based 
spintronics, supposing that such a technology will be fea- 
sible at all in carbon nanosy stems. 

We hope that our results might also lead to an incresing 
interest in boron based nanomaterials, whose polymor- 
phism is at least as rich as the polymorphism of carbon 
nanomaterials^^^^. Our present results indicate that in- 
teresting functionalities might actually be gained from 
the combination of carbon and boron nanosystems, and 
this might be a promising new field to be explored in the 
future, besides pure graphene. 
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